NASA Technical Memorandum 110162 




Formulation of an Improved Smeared 
Stiffener Theory for Buckling Analysis 
of Grid-Stiffened Composite Panels 


Navin Jaunky and Norman F. Knight, Jr. 
Old Dominion University , Norfolk, Virginia 

Damodar R. Ambur 

Langley Research Center, Hampton, Virginia 


(NASA-TM-110162) FORMULATION OF AN N95-30341 

IMPROVED SMEARED STIFFENER THEORY 
FOR 6UCKLING ANALYSIS OF 

GRID-STIFFENED COMPOSITE PANELS Unclas 

(NASA. Langley Research Center) 

16 p 

G 3/24 00 5 5 5 68 


June 1995 


National Aeronautics and 
Space Administration 
Langley Research Center 
Hampton, Virginia 23681-0001 



FORMULATION OF AN IMPROVED SMEARED 

STIFFENER THEORY FOR BUCKLING ANALYSIS 
OF GRID-STIFFENED COMPOSITE PANELS 


Navin Jaunky and Norman F. Knight, Jr. 
Old Dominion University 
Norfolk, VA 23529-0247 

Damodar R. Ambur 
NASA Langley Research Center 
Hampton, VA 23681-0001 


ABSTRACT 

A smeared stiffener theory for stiffened panels is presented that includes skin- 
stiffener interaction effects. The neutral surface profile of the skin-stiffener combi- 
nation is developed analytically using the minimum potential energy principle and 
statics conditions. The skin-stiffener interaction is accounted for by computing the 
stiffness due to the stiffener and the skin in the skin-stiffener region about the neu- 
tral axis at the stiffener. Buckling load results for axially stiffened, orthogrid, and 
general grid-stiffened panels are obtained using the smeared stiffness combined with 
a Rayleigh-Ritz method and are compared with results from detailed finite element 
analyses. 

INTRODUCTION 

In aircraft structures, structural efficiency dictates that most primary structures 
be of stiffened construction. The advent of high-performance composite materials 
combined with low-cost automated manufacturing using filament-winding and tow- 
placement techniques has made grid- stiffened structural concepts a promising al- 
ternative to more traditional stiffened structural concepts. Their damage tolerant 
characteristics (Ref. 1) and stiffness tailoring potential (Ref. 2) make grid-stiffened 
structures attractive for structural applications. 

An aircraft in flight is subjected to air loads which are imposed by maneuver and 
gust conditions. These forces on a structured panel that result from these external 
loads are shown in Figure 1. These internal loads, which depend on the location of 
the panel in an aircraft structure, may result in overall panel buckling, local buckling 
of the skin between stiffeners, and stiffener crippling. Hence, an efficient and accurate 
buckling analysis method for general grid- stiffened panels subjected to combined in- 
plane loading is needed in order to design grid-stiffened structural panels for different 
locations in fuselage and wing structures. 

Most of the research work on stiffened panels presented in the literature addresses 
axially stiffened panels subjected to compression. A limited amount of work has been 



reported on stiffened panels subjected to combined in-plane loading. Axially stiffened 
panels subjected to axial compression and in-plane shear was considered by Stroud, 
et al. (Ref. 3). Gendron and Gurdal (Ref. 4) considered grid-stiffened composite 
cylindrical shells subjected to axial compression and torsional shear. The modeling 
approaches generally used in the analysis of stiffened panels include the discrete ap- 
proach (Ref. 5), the branched plate and shell approach (Refs. 3 and 4), and the 
smeared stiffener approach (Refs. 6-9). In the discrete stiffener approach, stiffen- 
ers are modeled as lines of axial bending and torsional stiffnesses on the skin. This 
approach is difficult to use when the panel is stiffened in more than two directions 
and when the stiffener is not symmetric about the skin mid-surface. The branched 
plate and shell approach is more flexible and more accurate and usually involves the 
use of finite element analysis (Ref. 4). However, the detailed spatial discretization 
of the finite element model is tedious, and the solution is computationally expensive. 
In the smeared stiffener approach, the stiffened panel is converted mathematically to 
an unstiffened uniform thickness panel with equivalent orthotropic stiffnesses. These 
equivalent or smeared stiffnesses can be used in a Rayleigh-Ritz method to solve for 
buckling loads of the stiffened panel. The smeared stiffener approach is computa- 
tionlly efficient to execute and can easily account for stiffeners in any direction. The 
smeared stiffener approach is applicable in general to stiffened panels where the local 
buckling load is equal to or greater than the global buckling load. This approach 
for preliminary design is consistent with the aeronautical design philosophy where a 
buckling-resistant design is the design goal. 

In Refs. 8 and 9, a first-order shear- deformation theory (FSDT) has been used 
for developing an analysis tool based on a smeared stiffener approach. As observed 
in Refs. 3 and 9, the traditional or conventional smeared stiffener approach may 
overestimate the buckling load of stiffened panels in a certain range of geometric 
parameters because the traditional smeared stiffener approach does not account for 
local skin- stiffener interactions. This effect should be included in an improved smeared 
stiffener approach to make the approach a more reliable tool for the analysis and 
design of grid- stiffened panels. This paper describes an approach to incorporate the 
effects of local skin-stiffener interaction into a smeared stiffener theory and presents 
numerical results for panel buckling loads from the traditional and improved smeared 
stiffener theories. 

ANALYTICAL APPROACH 

The approximate stiffness added by a stiffener to the skin stiffness can be deter- 
mined by locating the position of the neutral surface in a skin- stiffener combination. 
The location of the neutral surface is determined theoretically through a study of 
the local stress distribution near the skin- stiffener interface similar to the approach 
presented in Ref. 10 for a panel with a blade stiffener. However, the study presented 
in Ref. 10 does not provide a general solution that is applicable to all classes of 
symmetric laminates. 
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A grid-stiffened panel may be considered to be an assembly of repetitive units 
or unit cells (see Figure 2). Any stiffener segment in the unit cell may be isolated 
in a semi-infinite skin- stiffener model as shown in Figure 2 for a diagonal stiffener. 
An approach for obtaining the stress distribution in a semi-infinite stiffened panel is 
outlined below. 

The average membrane stresses in the local coordinate system of the semi-infinite 
stiffened panel model are obtained by combining the constitutive relations with the 
strain compatibility equations and the use of a stress function approach. As a result , 
the following fourth-order partied differential equation is obtained. 
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where A*- is given by [A t j/t]~ l , the Aij are the extensional stiffness coefficients of the 
skin and t is the thickness of the skin. Dividing Equation (1) by A*! and transforming 
the y coordinate by 77 = eo y results in 

B'F A;, &F .., (2 A\i + A\ t ) PF n J ,AU d'F 9 4 F _ 
dx 4 ° A\ 1 dx z dr) 0 Ah dx 2 drj 2 0 Ahdxdr ) 3 dr) 4 


where eo = [AJj/A^] 1 ^ 4 - This equation is solved by assuming that stresses decay 
rapidly as the distance, y , away from the stiffener centerline becomes large, that the 
stresses are localized near the stiffener, and that a symmetric loading condition exists 
along the stiffener. The membrane stress function is assumed to be of the form 


F = Real ( e imk ( x + ire oy ) ) = Real{e xmk{x+ ' rrt) ) 


(3) 


where k = m = 1,2,3, ... and r is an unknown, and x and y axe local coordinates 
in the semi-infinite model. Substituting this stress function into the fourth-order dif- 
ferential equation results in a quartic equation in terms of the unknown r. The roots 
of the quartic equation are computed using subroutine CXPOLY from the Mathe- 
matical and Statistical Software (Ref. 11) at NASA Langley Research Center. The 
roots of the quartic equation occur as two pairs of complex numbers given by 

r= t±rm+irn (4) 

l ±t\R2 + It ii 

The membrane solution corresponds to the root with the largest magnitude of the 
real part for r and is developed as follows 

F lm = Real [ e * mfe l iB +*( r Ji+*>7 h) ] = e - rRmkr >cos[mk{x - rj 17)], for rj > 0 

F 2m = Real[e imk[x+i( -- TR+iTl)n] ] = e- rRTnkr, cos[mk(x + riT})], for 77 < 0 
F m = ^(F lm + Fi m ) = A m e~ TnkecTR ( y - t ‘/ 2 }C os[mkx]C os[mkeori(y - tj 2)] (5) 

where tr and 77 are the real and imaginary parts of the root, respectively, t s is the 
thickness of the stiffener, and A m are the unknown coefficients to be determined. 
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A similar approach is taken for the bending solution using the fourth-order partial 
differential equation for the out-of-plane deflection in terms of the local coordinate 
system. That is, 


d^W_ . £>16 d 4 W 2 ( 2£>12 + £> 66 ) dSv 3 £>26 d A w d 4 w _ 

dx A b Dudx 3 dq b Du dx 2 dt] 2 eb Dudxdq 3 dq 4 


( 6 ) 


where £>;j are the bending stiffness coefficients of the skin, e b = [£>n/£> 22 ] 1//4 and 
q = e b y. The solution for the out-of- plane deflection is obtained by assuming that the 
out-of-plane deflection decays as y becomes large and that the loading is symmetric 
along the stiffener. The out-of-plane deflection is assumed to be of the form 


w = Real{e' mk(x+ir '^) 


( 7 ) 


which on substitution into Equation (6) gives another quartic equation in r. The 
solution for the out-of-plane displacement corresponds to the root with the smallest 
non-zero magnitude of the real part for r and is developed as follows 

wi m = for 7] >0 

w 2m = f 0 r jj < 0 

-j- C m W2m ~~t"" C m W2m\ 

W m = e -mfeekrjlt ( v-t ' /2 ){£? m sin[mke b r Ib (y - t,/ 2)] 

-\-Cm cos[mke b ri b (y - tj 2)]} Cos[mkx ] (8) 


where rja, and rj b are the real and imaginary parts of the root, respectively, and B m 
and C m are the unknown coefficients to be determined. 

These two solutions (Equations (5) and (8) ) are valid near the skin-stiffener 
interface but not within the stiffener itself (i.e., y > t s j 2). It is assumed that, since 
the stiffener is thin, the strain within the stiffener is approximately equal to the 
strain at the edge of the stiffener (at y = t,/2). The total strain energy, Ur-, of 
the skin-stiffener combination is developed next from expressions for the out-of-plane 
deflection, w m , and the membrane stress function, F m . The total strain energy is 
obtained by evaluating the following integrals 


1. The strain energy of the skin is 


Utkin = f f ( {eo} T [>M{eo} + {«} T [Aj]{«} ) dxdy (9) 

Jt,/2 J —L 

where {e 0 } = 7°} are the membrane strains and {«} = {k x k v n xy } are 

the curvatures. 

2. The strain energy of the stiffener is 


1 ft,/2 />— 1/2 fL 

U stiff = r / Qu(e° x +zk x ) 2 y=ti/2 dxdzdy ( 10 ) 

2t Jtg/2 J — (t/2+h) J— L 

where Qu is the longitudinal modulus of the stiffener. 
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3. The strain energy of the skin attached to the stiffener is 

Uattach = "f J_ L ( ^ll( £ x)y=t./2 + )^ x (H) 

Hence, the expression for the total strain energy, Ut, is obtained by summing these 
contributions that results in 

U T = C A A 2 m + C b BI + C c C 2 m + C AC A m C m + CcsCmBm (12) 

where the coefficients C A , Cb, Cc, C A c and Ccb are obtained by evaluating the strain 
energy integrals. 

The total bending moment developed at any cross section perpendicular to the 
longitudinal axis of the stiffener for the symmetric case can be represented by the 
series 

OO 

M = ^ M^(y) cos(mkx) (13) 

m=l 

From statics, the normal stresses over the cross section of plate-stiffener combination 
must satisfy the following conditions 

fit / 2 too fit / 2 

2/ / cr x dzdy + t , Q n ( e° x + zk x ) y=t ,/ 2 dz - 0 (14) 

Jt/2 h.j 2 J-{t/2+h) 

fit / 2 fi OO fit/2 

2 / za x dzdy + t s zQu( + zk x ) y =tJ 2 dz 

Jt/2 Jt.l 2 7-(t/2+/i) 

= ^2 M m cos(mkx ) (15) 

771 = 1 

where t is the total thickness of the skin, t , is the total thickness of the stiffner, 
h is the height of the stiffener above the outer surface of the skin, and <r x is the 
normal stress distribution over the cross-section. Evaluating the integrals defined by 
Equations (14) and (15) results in the following relations after neglecting coefficients 
of sin(mkx) which are due to the Aie and Di6 terms in the extensional and bending 
stiffness matrices, respectively. 


5'll-4.m + SlzCm — 0 

S 21 A m + S 22 B m + S 23 C m = M m (16) 

Using Equations (16), the following expressions for B m and C m are obtained in terms 
of A m and M m 


Cm = S* u A m 

B m ■ A m + S2 2 M m 


(17) 
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where 


S* n = -5 n /S 13 , S* 21 = (5 n 5 23 - 5 21 5 13 ) / S 13 S 22 

and S* 2 = 1/S 22 (18) 

Equations (17) are substituted into Equation (12) which is minimized with respect 
to A m to yield 

Am = -V M M m /V A (19) 

where 


v A = 2 (c a + c b (s ; 1 ) 2 + c c (s; 1 ) 2 CcbSuS ; 1 

v M = 2 C b s; 1 s ; 2 + CcbS- 11 s; 2 (20) 


Using Equation (19), B m and C m can be expressed in terms of M m , Vm and V A , with 
as the only unknown. 


C - 5* a - -S* V mM - 

VA 

VmM„ 

■^21 \T 


B m — A m "I - S 22 M m — 


"t" S 22 M m 


( 21 ) 


The expression for axial strain in the skin-stiffener combination is obtained from the 
stress function, F , and the out-of-plane deflection, w, which is given by 

1 d 2 F d 2 w 
£ * = [An It) Sy> ~ 

Substituting for A m , B m and C m in Equation (22) from Equations (19) and (21) and 
on solving for the value of z for which e x is zero, an expression for the neutral surface, 
Z'(y) is obtained. Only one term (m = 1) in the series expansion is used to obtain 
the expression for Z'(y). 
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The expression for Z\y) is independent of M m and the axial distance x . Since the 
expression for Z r involves Exp[mk(ei,r} jy> — e 0 rji)(y — f 2 /2)], the choice of roots for 
the solution of the stress function, jF, and the out-of-plane deflection, tu, ensures that 
the neutral surface Z\y ) decays as the distance away from the centerline, y, becomes 
large. Finally, the shift in the neutral surface at the stiffener is obtained by setting 
y = t a / 2 in the expression for Z\y). 

7 4(r 2 R-r}) 

(A 11 /t)S* 11 

A typical profile of the neutral surface for a skin-stiffener combination is shown in 
Figure 3. The distance y * represents the distance from the centerline of the stiffener 
to the point where the neutral surface coincides with the mid-surface of the skin. The 
average distance of the neutral surface over the distance y* is Z*. The quantities y* 
and Z* are obtained numerically. The correction to the smeared stiffnesses due to 
the skin- stiffener interaction is introduced by computing the stiffness of the stiffener 
and the skin segment directly contiguous to it according to the following criteria. 

1. If y* < f/4, then the reference surface for the stiffener is Z„. 

2. If y* > f/4, then the reference surface for the stiffener is Z*. 

In either case, the reference surface of the skin is taken to be its mid-surface. 

NUMERICAL RESULTS 

Three stiffened panels with different stiffener configurations and simply-supported 
boundary conditions are used as examples for the present analytical approach. Panel 
1 is an axially-stiffened panel, Panel 2 is an orthogrid-stiffened panel, and Panel 3 is 
an example for a general grid- stiffened panel. Finite element analyses of these three 
panels have been conducted to verify the results for the present analytical approach. 
The finite element analysis codes STAGS (Ref. 12) and DIAL (Ref. 13) have been 
used for this purpose. In the STAGS finite element model, a nine-node shear-flexible 
element (i.e., STAGS element 480) is used while an eight-node isoparametric shear 
flexible element is used in the DIAL model. Finite element analysis results for all 
panels indicate that the panels buckle globally under the applied in-plane loading 
conditions. 

Panel 1 

Panel 1 is 30.0-in. (762-mm.) long and 30.0-in. (762-mm.) wide with axial 
stiffeners only. The stiffener height and thickness are 1.86958 in. (47.5 mm.) and 
0.20084 in. (5.1 mm.), respectively. The unit cell is 30.0-in. (762-mm.) long and 
10.0-in. (25.4-mm.) wide (see Figure 4). The skin ply stacking sequence is [45/ — 
45/ — 45/45/0/90], with thicknesses of 0.00637 in. (0.16 mm.) for the 45° and —45° 
plies, 0.0249 in. (0.63 mm.) for the 0° plies and 0.0416 in. (1.05 mm.) for the 90° 
plies. The stiffener ply stacking sequence is [45/ — 45/ — 45/45/0], with thicknesses 
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of 0.00823 in. (0.21 mm.) for the 45° and —45° plies and 0.0675 in. (1.71 mm.) for 
the 0° plies. The nominal ply mechanical properties used are: longitudinal modulus 
= 19.0 Msi (131.16E03 MPa); transverse modulus = 1.89 Msi (13.04E03 MPa); shear 
modulus = 0.93 Msi (6.42E03 MPa) and major Poisson’s ratio = 0.38. 

The four panel load cases considered are shown in Table 1. The STAGS analysis 
results are compared with solutions from the smeared stiffener approach without skin- 
stiffener interaction effects included (the traditional approach) and with skin-stiffener 
interaction effects included (the present approach). It can be seen that the value of 
Z n for the axial stiffener is not small compared to the height of the stiffener. The 
result obtained from the traditional approach is in good agreement with the STAGS 
analysis result for the case of axial compression and the result from present approach 
is less than the STAGS analysis result by 7.5 percent. For the other load cases shown 
in the Table, the results obtained by the traditional approach are greater than those 
of STAGS by 8 to 13 percent and those of the present approach are in good agreement 
with the STAGS results. 

Panel 2 

Panel 2 is 60.0-in. (1524-mm.) long and 36.0-in. (914.4-mm.) wide with axial 
and transverse stiffeners only. The stiffener height and thickness are 0.5 in. (12.7 
mm ) and 0.12 in. (3.0 mm.), respectively. The unit cell is 20.0-in. (508-mm.) 
long and 9.0-in. (228.6-mm.) wide (see Fig. 5). The skin ply stacking sequence 
is [45/ — 45/90/0] s and each ply thickness is 0.008 in. (0.20 mm.). The stiffener is 
made of material with 0° orientation. The nominal ply mechanical properties used 
are: longitudinal modulus = 24.5 Msi (169.13E03 MPa); transverse modulus = 1.64 
Msi (11.32E03 MPa); shear modulus = 0.87 Msi (6.0E03 MPa) and major Poisson’s 
ratio = 0.3. 

The panel buckling response when subjected to four loading conditions is indi- 
cated in Table 2. The DIAL analysis results are compared in Table 2 with solutions 
from the smeared stiffener approach without skin-stiffener interaction effects and with 
skin- stiffener interaction effects. The value of Z n for the transverse stiffener is not 
small compared to the height of the stiffener. The results obtained using the tradi- 
tional approach overestimate the DIAL analysis result by 12.6 percent for the axial 
compression load case, by 4.0 percent for the transverse compression load case, and 
by 8.4 percent for the combined load cases. Results from the present approach agree 
with the DIAL analysis results except for the transverse compression load case where 
the present result is 5.2 percent less than the DIAL analysis result. 

Panel 3 

Panel 3 is 56.0-in. (1422.4-mm.) long and 20.0-in. (508-mm.) wide with trans- 
verse and diagonal stiffeners only. The stiffener height and thickness are 0.276 in. 
(7.0 mm.) and 0.1125 in. (2.86 mm.), respectively. The unit cell dimensions for this 
panel are 7.0 in. (177.8 mm.) in length and 5.0 in. (127 mm.) in width (see Fig. 
6). The skin stacking sequence is [45/90/ - 45] „ and each ply thickness is 0.008 in. 
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(0.20 mm.). The stiffener for this case is also made of 0° material. The nominal ply 
mechanical properties used are: longitudinal modulus = 24.5 Msi (169.13E03 MPa); 
transverse modulus = 1.64 Msi (11.32E03 MPa); shear modulus = 0.87 (6.0E03 MPa) 
Msi and major Poisson’s ratio = 0.3. 

The panel was analyzed for the three load conditions shown in Table 3. The 
DIAL analysis results are compared with results from the smeared stiffener approach 
without skin-stiffener interaction effects and with skin-stiffener interaction effects in 
Table 3. For this panel, the values of Z n are small compared to the height of the 
stiffener. The results obtained from the traditional approach are approximately 11 
percent greater than the DIAL analysis results, and the results obtained using the 
present approach are approximately 6.5 less than the DIAL analysis results. For 
this panel, the results obtained using the present approach are conservative since the 
contribution of stiffness terms A\e and Die in the expression for a x are not small and 
influence the neutral surface profile position for the diagonal stiffener. 

CONCLUDING REMARKS 

An improved smeared stiffener theory that includes skin-stiffener interaction ef- 
fects has been developed. The skin-stiffener interaction effects are introduced by 
computing the stiffness of the stiffener and the skin at the stiffener region about the 
neutral axis at the stiffener. The neutral surface profile for the skin-stiffener combi- 
nation is obtained analyticaly through a study of the local stress distribution near 
the skin-stiffener interface. 

The results from the numerical examples considered suggest that skin-stiffener 
interaction effects should be included in the smeared stiffener theory to obtain good 
corelation with results from detailed finite element analyses. In a few cases the present 
analysis appears to underestimate the buckling load by 5 to 7 percent. In spite of this 
limitation, the smeared stiffener theory with skin-stiffener interaction effects included 
is still a useful preliminary design tool and results in buckling loads that are more 
accurate than the results from the traditional smeared stiffener approach. 
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Table 1: Results for axially stiffened panel (Panel 1). 


X-stiffener: Z n = -0.4386 in., Z* = -0.1020 in., y* = 4.7512 in. 


Z n = -11.14 mm., Z* = -2.59 mm., y* = 120.67 mm. 




Critical Eigenvalue 

N x 

lbs/in. 

N xy STAGS 

lbs/in. 

'< Traditional 
Approach 

Present 

Approach 

1000 

(175.34) 

0 9.9636 

(0) 

9.9659 

9.2135 

0 

(0) 

1000 6.3016 

(175.34) 

6.7985 

6.3483 

1000 

(175.34) 

1000 4.9512 

(175.34) 

5.6018 

4.9491 

500 

(87.67) 

1000 5.5023 

(175.34) 

6.2007 

5.5838 

Numbers within parentheses indicate loading in N/mm. 


Table 2: Results for orthogrid panel (Panel 2). 

X-stiffener: Z n = -0.0949 
Z n = -2.41 
Y-stiffener: Z n = -0.1295 
Z n = -3.29 

in., Z* = -0.0165 in., 
mm., Z* — -0.42 mm. 
in., Z* = -0.0177 in., 
mm., Z* = -0.45 mm. 

y* = 0.0280 in. 

, y* - 0.71 mm. 
y * = 0.0131 in. 

, y* = 0.33 mm. 



Critical Eigenvalue 

N x 

lbs/in. 

Ny N xy 

lbs/in. lbs/in. 

DIAL Traditional 
Approach 

Present 

Approach 

400 

(70.14) 

0 0 

(0) (0) 

0.7909 0.8903 

0.8161 

0 

(0) 

200 0 

(35.07) (0) 

0.6281 0.6536 

0.5956 

400 

(70.14) 

200 0 

(35.14) (0) 

0.3504 0.3799 

0.3463 

400 

(70.14) 

200 50 

(35.14) (8.77) 

0.3500 0.3796 

0.3458 


Numbers within parentheses indicate loading in N/mm. 
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Table 3: Results for grid-stiffened panel (Panel 3). 

Y-stiffener: Z n = -0.0135 in., Z* = -0.0043 in., y * = 2.3636 in. 

Z n = -0.34 mm., Z* = -0.11 mm., y* - 60.0 mm. 
D-stiffener: Z n = -0.0698, Z * = -0.0349 in., y* = 0.0239 in. 




= -1.77 

mm., Z* 

-0.89 mm., y* 

= 0.61 mm. 





Critical Eigenvalue 

N x 

N y 

N xy 

DIAL 

Traditional 

Present 

lbs/in. 

lbs/in. 

lbs/in. 


Approach 

Approach 

o o 
o 

400 

(70.14) 

0.0 

(0.0) 

0.3290 

0.3646 

0.3045 

0.0 

(0.0) 

400 

(70.14) 

300 

(52.60) 

0.3224 

0.3595 

0.3008 

100 

(17.53) 

400 

(70.14) 

300 

(52.60) 

0.3121 

0.3486 

0.2917 


Numbers witbin parentheses indicate loading in N/mm. 
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TRANSPORT 



Figure 1. Aircraft structural applications showing internal forces. 



stiffened 



Figure 2. Semi-infinite plate model for a skin-stiffener element. 



Figure 3. Typical profile of neutral surface for a 
skin-stiffener element. 
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Figure 4 . Axially stiffened panel 



D-stiffener 
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